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ABSTRACT 

In the solar system, quasi-satellites move in a 1:1 mean motion resonance going around their 
host body like a retrograde satellite but their mutual separation is well beyond the Hill radius 
and the trajectory is not closed as they orbit the Sun not the host body. Although they share 
the semi-major axis and the mean longitude of their host body, their eccentricity and incli- 
nation may be very different. So far, minor bodies temporarily trapped in the quasi-satellite 
dynamical state have been identified around Venus, Earth, the dwarf planet (1) Ceres, the large 
asteroid (4) Vesta, Jupiter, Saturn and Neptune. Using computer simulations, Tiscareno and 
Malhotra have predicted the existence of a small but significant population of minor bodies 
moving in a 1:1 mean motion resonance with Pluto. Here we show using TV-body calculations 
that the Plutino 15810 (1994 JRi) is currently an accidental quasi-satellite of Pluto and it will 
remain as such for nearly 350,000 yr. By accidental we mean that the quasi-satellite phase is 
triggered (or terminated) not by a direct gravitational influence in the form of a discrete close 
encounter but as a result of a resonance. The relative mean longitude of the Plutino circulates 
with a superimposed libration resulting from the oscillation of the orbital period induced by 
the 2:3 mean motion resonance with Neptune. These quasi-satellite episodes are recurrent 
with a periodicity of nearly 2 Myr. This makes 15810 the first minor body moving in a 1:1 
mean motion resonance with Pluto and the first quasi-satellite found in the trans-Neptunian 
region. It also makes Pluto the second dwarf planet, besides Ceres, to host a quasi-satellite. 
Our finding confirms that the quasi-satellite resonant phase is not restricted to small bodies or- 
biting major planets but is possible for dwarf planets/asteroids too. Moreover, 15810 could be 
considered as a possible secondary target for NASA's Pluto-Kuiper Belt Mission New Hori- 
zons after the main Pluto flyby in 2015. This opens the possibility of studying first hand and 
for the first time a minor body in the quasi-satellite dynamical state. 

Key words: celestial mechanics - planets and satellites: individual: Pluto - asteroids: indi- 
vidual: 15810 (1994 JRI) - Solar System: general- minor planets, asteroids 



1 INTRODUCTION 

Quasi-satellites are minor bodies that share the semi-major axis and 
the mean longitude of their host body but may have different ec- 
centricity and inclination. They move like a retrograde satellite but 
their separation is well outside the Hill sphere of the host body and 
the trajectory is not closed; therefore, they are not bound satellites 
as they orbit the Sun not the host body. The theory behind these 
remarkable objects was first studied in 1913 (Jackson 1913) but 
the topic was largely neglected by the scientific community until 
the end of the XXth century when co-orbital motion of planets and 
asteroids started receiving more attention on theoretical grounds 
(Mikkola & Innanen 1997; Wiegert, Innanen & Mikkola 2000). So 
far, such quasi-satellites have been found around Venus (Mikkola et 
al. 2004), Earth (Wiegert, Innanen & Mikkola 1997; Connors et al. 
2002; Connors et al. 2004; Brasser et al. 2004; Christou & Asher 
2011), the dwarf planet (1) Ceres and the large asteroid (4) Vesta 
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(Christou & Wiegert 2012), Jupiter (Kinoshita & Nakai 2007; Wa- 
jer & Krolikowska 2012), Saturn (Gallardo 2006) and Neptune (de 
la Fuente Marcos & de la Fuente Marcos 2012a). 



If quasi-satellites have been discovered orbiting rocky bodies 
in the inner solar system and theory predicts long-term stability 
for quasi-satellite orbits in the outer solar system, one may wonder 
whether a quasi-satellite could orbit an object like the dwarf planet 
Pluto. The existence of a small but significant population of mi- 
nor bodies experiencing co-orbital resonant behaviour with respect 
to Pluto in the form of libration or slow circulation of the relative 
mean longitude has been predicted in the context of chaotic diffu- 
sion of trans-Neptunian objects (Yu & Tremaine 1999; Tiscareno & 
Malhotra 2009). Such objects may experience relatively close, low 
velocity encounters with Pluto which translate into comparatively 
large perturbational effects if they occupy the Kozai resonance (Tis- 
careno & Malhotra 2009). Using data from the JPL's HORIZONS 
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systerrQ we performed a numerical survey looking for minor bodies 
relatively close to Pluto during the next few years in order to find 
candidates for possible co-orbital resonant behaviour. Our search 
resulted in a promising candidate, 15810 (1994 JRi). We identi- 
fied 15810 as currently located at about 3.1 AU from Pluto and 
slowly approaching to its peripluto at 2.7 AU within a timeframe of 
5 years. The object was originally discovered with the 2.5 m Isaac 
Newton Telescope at La Palma on May 12, 1994 (Irwin et al. 1994; 
Irwin, Tremaine & Zytkow 1995) and it has a diameter of 251 km 
(Irwin et al. 1995). Its orbit, which is quite reliable, has been com- 
puted using 43 observations with an arc- length of 2236 day!0 and 
its UBVRI colors have also been obtained (Barucci et al. 1999). 

In this Letter and with the help of TV-body calculations, we 
show that the Plutino asteroid 15810 (1994 JRi) currently follows 
a quasi-satellite path around Pluto. This Letter is organized as fol- 
lows: in Section 2, we briefly outline our numerical model. Section 
3 presents and discusses our results. Our conclusions are summa- 
rized in Section 4. 



2 NUMERICAL MODEL 

In order to investigate the orbital evolution of 15810 and its pos- 
sible co-orbital resonant behaviour with respect to Pluto we have 
performed TV-body calculations in both directions of time. The nu- 
merical integrations of the orbit of Plutino 15810 presented here 
were computed with the Hermite integrator (Makino 1991; Aarseth 
2003), in a model solar system which included the perturbations by 
the eight major planets (Mercury to Neptune) and treat the Earth 
and the Moon as two separate objects, it also includes the barycen- 
tre of the dwarf planet Pluto-Charon system and the three largest 
asteroids, (1) Ceres, (2) Pallas and (4) Vesta. The standard version 
of this direct TV-body code is publicly available from the IoA web 
sitfi additional details can be found in de la Fuente Marcos & de la 
Fuente Marcos (2012b). Results in the figures have been obtained 
using initial conditions (positions and velocities in the barycentre 
of the solar system referred to the JD2456200.5 epoch) provided by 
the JPL Horizons system (Giorgini et al. 1996; Standish 1998). In 
addition to the calculations completed using the nominal orbital el- 
ements in TableQ]we have performed 100 control simulations using 
sets of orbital elements obtained from the nominal ones using the 
accepted uncertainties (3er). The sources of the Heliocentric Keple- 
rian osculating orbital elements of 15810 are the JPL Small-Body 
DatabastiJ and the AstDyS-2 portafl 



3 ORBITAL EVOLUTION 

The heliocentric orbits of both 15810 and the barycentre of the 
Pluto-Charon system (panels A and B) as well as the quasi-satellite 
motion features (panels C and D) are shown in Fig. Q] The motion 
of 15810 from 2012 to 17,012 (the origin of time = JD2456200.5, 
2012-Sep-30.0) shows quasi-satellite loops ("corkscrew" orbits) as 
viewed from above Pluto (C) and from a point outside Pluto's orbit 
looking past Pluto in towards the Sun (D), in a frame of reference 

1 http://ssd.jpl.nasa.gov/7horizons 

2 http://www.minoiplanetcenter.net/db_search/show_object ?object_id= 
1994+JR1 

3 http://www.ast.cam.ac.uk/~sverre/web/pages/nbody.htm 

4 http://ssd.jpl.nasa.gov/sbdb.cgi 

5 http://hamilton.dm.unipi.it/astdys/ 



Table 1. Heliocentric Keplerian orbital elements of 15810 (1994 JRi) used 
in this research. Values include the 1-cr uncertainty. (Epoch = JD2456200.5, 
2012-Sep-30.0; J2000.0 ecliptic and equinox. Source: JPL Small-Body 
Database and AstDyS-2 portal.) 



Semi-major axis, a 


= 39.24±0.02 AU 


Eccentricity, e 


= 0.1143±0.0003 


Inclination, i 


= 3.8032±0.0002 ° 


Longitude of ascending node, Q 


= 144.753±0.011 ° 


Argument of perihelion, ui 


102.1±0.2° 


Mean anomaly, M 


= 24.42±0.12° 



revolving with Pluto. Each loop takes one Plutonian year, 247.7 yr. 
In order to further study the resonant properties of the path shown in 
Fig. H] C-D panels, let us define the relative deviation of the semi- 
major axis from that of Pluto by a = (a — ap)/ap, where a and 
ap are the semi-major axes of 15810 and Pluto respectively; and 
also the relative mean longitude A — Ap, where A and Ap are the 
mean longitudes of 15810 and Pluto respectively. In the top panel, 
Fig. [2] the evolution of a as a function of A — Ap during the time 
interval (-25,000, 100,000) yr is displayed. The short period fluctu- 
ations are associated to the period of Pluto. In principle, the secular 
motion is a quasi-harmonic oscillation of the variables A — Ap and 
a; this is the main feature of the quasi-satellite motion (Mikkola et 
al. 2006). In the middle panel, the mean longitude of 15810 relative 
to Pluto is displayed. The object currently librates asymmetrically 
around 0° with amplitude 40° -50° and a period of about 20,000 yr 
that coincides with the libration period of the main resonant angle 
of a 2:3 mean motion resonance with Neptune (see below). Here 
by amplitude we mean the difference between the maximum and 
the minimum values of the relative mean longitude in a period. The 
object slowly drifted from the neighbourhood of L4 into the quasi- 
satellite dynamical state and, in the future, it will move towards L5 . 
L4 is the Lagrange point 60° ahead of Pluto, L5 is the Lagrange 
point trailing Pluto 60°; in general, motion around the Lagrange 
triangular points follows a tadpole orbit and the objects in these 
orbits are called Trojans. 

Our object, however, does not follow a classical tadpole or 
horseshoe behaviour after leaving its quasi-satellite path due to the 
slow circulation likely induced by the close approaches with Pluto. 
The actual quasi-satellite phase lasts nearly 350,000 yr. Before and 
after the actual quasi-satellite state, the object does not follow the 
classical compound orbits resonating between the Trojan and quasi- 
satellite dynamical states typical of other quasi-satellites (e.g. 2002 
VEs8, Mikkola et al. 2004). These compound orbits have been de- 
scribed on theoretical grounds (Namouni 1999; Namouni, Christou 
& Murray 1999). The evolution of the semi-major axis over the 
plotted periods (Figs.[2]and[3} remains fairly stable and its average 
value (proper semi-major axis), 39.4518 AU, is very close to that of 
Pluto, 39.4477 AU. All our simulations suggest that the object has 
remained in the quasi-satellite phase for nearly 100,000 yr. In the 
future, 15810 will leave the quasi-satellite path slowly drifting to- 
wards the L5 Lagrangian point. This will happen about 250,000 yr 
from now. The long-term evolution of the orbital elements of 15810 
displayed in Fig. [3] shows that the eccentricity exhibits ~0.5 Myr 
periodic variations, likely the result of an unidentified secular reso- 
nance. The inclination also oscillates but the argument of perihelion 
does not librate about a value of 90° like in the case of Pluto, it cir- 
culates with a period of ~0.7 Myr. Therefore, 15810 is not submit- 
ted to a Kozai's secular resonance (Kozai 1962). Close encounters 
with Pluto have a periodicity of ~2 Myr. The difference between 
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Figure 1. Orbital evolution of Plutino 15810 (1994 JRi) in the time interval 
(0, 15,000) yr. A-B. The orbits of 15810 (thick line) and Pluto (the barycen- 
tre of the Pluto-Charon system) when seen in the heliocentric J2000 ecliptic 
frame of reference with the x-axis aligned towards the vernal equinox. In A 
they are seen projected from the direction of the north ecliptic pole (ecliptic 
plane). In B the orbits are seen projected from the direction of the vernal 
equinox. The relatively large difference between the orbital inclination of 
15810 (3.80°) and Pluto (17.14°) is evident in the figure. C-D. The orbit of 
15810 in plutocentric coordinates co-rotating with Pluto (the barycentre of 
the Pluto-Charon system). The path oscillates in such a way that when the 
relative mean longitude librates around 0° (see Fig.0, Pluto remains inside 
the path of 15810. The orientation of the asteroid's orbit allows the path to 
overlap the position of Pluto without any danger of collision although close 
approaches are certainly possible (see the text). In all the figures the origin 
of time = JD2456200.5, 2012-Sep-30.0. 



158 10's and Pluto's mean longitudes also circulates every -~2 Myr. 
The similar periodicity strongly suggests that the close approaches 
with Pluto are responsible for the circulation of the relative mean 
longitude. If we repeat the calculations assuming a negligible mass 
for the Pluto-Charon system, the librating component of the relative 
mean longitude (Fig. [2] middle panel) vanishes and just the circu- 
lating behaviour (but with longer period) remains: Pluto is clearly 
responsible for the observed quasi-satellite behaviour. 

The distance of 15810 from Neptune and Pluto is shown in 
Fig. E] The distance of 15810 from Neptune remains larger than 
1 1.5 AU (the distance from Uranus > 13.5 AU). This indicates that 
encounters with Neptune (or Uranus) do not cause the object to de- 
part from the quasi-satellite orbit in either direction of time. This 
is to be expected as 15810 is itself a Plutino, the object is in a 2:3 
mean motion resonance with Neptune: in the time period 15810 
completes two orbits, Neptune goes around the Sun three times. 
The resonance argument <f> = 3A — 2 A at — w, with vj — $7 + uj 
and Ajv the mean longitude of Neptune, librates about 180° with 
an amplitude of 90° (not 82° like in the case of Pluto), the libration 
period is ~20,000 yr. This resonance prevents close encounters be- 
tween 15810 and Neptune and also between Pluto and Neptune. 
The slow drift in relative mean longitude observed in Fig.|2]is not 
induced by Neptune but by Pluto; the object experiences repeated 
relatively close encounters with Pluto (more properly, the barycen- 
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Figure 2. Resonant evolution of Plutino 15810 (1994 JRi). (top) The rel- 
ative deviation of its semi-major axis from that of Pluto (the barycentre 
of the Pluto-Charon system), a, as a function of the difference between 
its mean longitude from that of Pluto (the barycentre of the Pluto-Charon 
system) during the time interval (-25,000, 100,000). The relative mean lon- 
gitude librates around 0° which is the signpost of the quasi-satellite be- 
haviour (Mikkola et al. 2006). (middle) Mean longitude relative to Pluto, 
A - A P , over the time interval (-300,000, 300,000). Plutino 15810 is cur- 
rently following a quasi-satellite path with the relative mean longitude li- 
brating asymmetrically around 0° . (bottom) Semi-major axis evolution. The 
~20,000 yr periodic variations induced by the 2:3 mean motion resonance 
with Neptune are observed in the evolution of the semi-major axis. Initial 
conditions (nominal orbit) are given in Table[T] 

tre of the Pluto-Charon system as in all this discussion) as seen in 
Fig. E] although initially the closest distances are well outside the 
Hill sphere of Pluto that has a radius of 0.0385 AU. Approximately 
125,000 yr into the past, 15810 had a number of periodic (10,000 
yr) periplutos in the range 1.6-1.8 AU. But in the future, in about 
13,000 yr, much closer approaches will be possible at 0.07 AU that 
is less than twice its Hill radius. The object will slowly drift into 
the area around the Ls Lagrange point. 

The object identified in this Letter is not the classical quasi- 
satellite found around objects moving in not very eccentric orbits 
like Venus. The high inclination and eccentricity characteristic of 
the orbit of Pluto add new complexity to an already challenging 
dynamical situation. In principle, it may be argued that 15810 is 
not even co-orbital with Pluto, making the analysis carried out here 
dynamically unjustified. Following Namouni (1999), the co-orbital 
region in the case of a host object is defined as \a — ao ^ tho, 
where tho is the radius of the Hill sphere of the host object and 
ao, its semi-major axis. In the case of the first bona fide quasi- 
satellite, 2002 VE 6 8, and for the JD2456200.5 epoch, the semi- 
major axes of the quasi-satellite and Venus are 0.7237 AU and 
0.7233 AU respectively, and the Hill radius of Venus is 0.0067 AU. 
Therefore, Namouni's criterion obviously indicates that 2002 VEsg 
is co-orbital with Venus. The same naive application of the crite- 
rion to 15810 gives 39.24 - 39.36 > 0.0385,4(7 in violation 
of the co-orbitality criterion. Therefore and using this approach, 
15810 is not even co-orbital with Pluto and the quasi-satellite mo- 
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Figure 3. Time evolution of various orbital elements. The orbital elements a 
(top panel), e (second to top), i (third to top) and the argument of perihelion 
(bottom panel) are displayed. Plutino 15810 (1994 JRi)'s orbital period 
oscillates about a mean value which is exactly 3/2 that of Neptune but its 
argument of perihelion circulates, it does not librate about a value of 90° 
like in the case of Pluto. Therefore, 15810 is not submitted to a Kozai's 
secular resonance (Kozai 1962). 



tion pointed out above is not more than a happy coincidence. How- 
ever, if we use the proper orbital elements instead of osculating Ke- 
plerian orbital elements at a particular epoch, the criterion becomes 
|39.4518 - 39.4477| < 0.0401^117. The proper elements clearly 
confirm the co-orbital nature (with Pluto) of 15810; Plutino 15810 
is truly, albeit somewhat accidental, a quasi-satellite of Pluto. By 
accidental we mean that the quasi-satellite phase is not triggered 
or terminated by direct gravitational interaction in the form of a 
discrete close encounter (with a certain planet like the Earth in the 
case of 2002 VEes) but as a result of a resonance. The relative 
mean longitude of the Plutino circulates with a superimposed libra- 
tion resulting from the oscillation of the orbital period induced by 
the 2:3 mean motion resonance with Neptune. These quasi-satellite 
episodes are recurrent with a periodicity of nearly 2 Myr. The be- 
haviour found in our calculations have been previously described 
in numerical simulations (Yu & Tremaine 1999; Tiscareno & Mal- 
hotra 2009) which predict that only 7% of objects is expected to ex- 
perience persistent circulation of the relative mean longitude with 
respect to Pluto; therefore, the identification of one of these unusual 
objects provides a useful constraint for models studying the dynam- 
ics of the outer solar system as well as giant planet migration. Re- 
garding the origin of 15810, its currently very stable orbit suggests 
that it is not relatively recent debris originated in collisions within 
Pluto's system but perhaps a primordial Plutino formed around the 
same epoch Pluto came into existence. The object studied here may 
be part of an outer solar system analogue to the population of Main 
Belt asteroids recently found co-orbiting with the dwarf planet (1) 
Ceres and the large asteroid (4) Vesta (Christou & Wiegert 2012). 
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Figure 4. The distance of Plutino 15810 from Neptune (top panel) and from 
Pluto (the barycentre of the Pluto-Charon system, bottom panel). The dis- 
tance of 15810 from Neptune remains larger than 11.5 AU as it moves 
in a 2:3 mean motion resonance with Neptune but the object experiences 
repeated encounters with Pluto (its barycentre). The Hill radius of Pluto, 
0.0385 AU, is also indicated. 

Although the figures have been computed using the nominal orbit 
in TableQ] the other simulations gave very similar results, over the 
time interval shown. 



4 CONCLUSIONS 

Pluto's system continues being a source of controversy, unan- 
swered questions and surprises more than 80 years after its discov- 
ery. Pluto's planethood demotion in August 2006 still stirs debate 
today and the recent finding of a fifth moon orbiting Pluto by the 
HSTjjjust confirms the unexpectedly complex nature of the system. 
In our work, we show that 15810 currently follows a quasi-satellite 
orbit relative to Pluto; therefore and besides having 5 regular satel- 
lites, Pluto has at least one quasi-satellite. This makes 15810 the 
first minor body found moving in a 1:1 mean motion resonance 
with Pluto and the first quasi-satellite found in the trans-Neptunian 
region of the solar system. It also makes Pluto the second dwarf 
planet, besides Ceres, to host a quasi-satellite. Our finding also 
confirms that the quasi-satellite resonant phase is not restricted to 
small bodies orbiting major planets but it is possible for dwarf plan- 
ets/asteroids too. We also provide a new and somewhat unexpected 
mechanism to land minor bodies into the quasi-satellite dynami- 
cal state. On the other hand, Plutino 15810 is a natural candidate 
for a spacecraft rendezvous mission in the framework of NASA's 
Pluto-Kuiper Belt Mission New Horizons that is going to complete 
a flyby with Pluto in 2015 and then continue to explore one or more 
nearby trans-Neptunian objects in the time frame 2016-2020. It is 



http://hubblesite.org/newscenter/archive/releases/solar- 
system/pluto/2012/32/ 
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moving in such an orbit that this object could be a good candi- 
date for a body dynamically related to the Pluto-Charon formation 
event, and the determination of the physical properties of its surface 
by spectroscopic observations could be interesting in that respect. 
If 15810 is selected for the extended mission it will open the pos- 
sibility of studying in detail for the first time a minor body in the 
quasi-satellite dynamical state. 
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